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Abstract 

00 In this study we report the effects of an external electromagnetic field on the collective properties 

^_^ of unmagnetized plasmas. The calculations are carried out in the semi-classical approximation, 

Q^ i.e., the electromagnetic field is treated classically and the electrons from a quantum mechanical 

^ viewpoint. The results show that the collective modes are damped away smoothly and in a smaller 

^ frequency range than those reported by previous studies. An exponential-like decay for the plasmon 



frequencies as a function of the external field amplitude is readily observed. The results of previous 



^ studies are successfully obtained . We also find that the single photon processes has a pronounced 

(N 

^^ effect on the decrease of the frequency range of modulation. 
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I. INTRODUCTION 

Quantum-mechanical tools have been used extensively to deal with classical plasma 
physics phenomena for a long time |lH3]. One of the main reasons for the success of such an 
approach is that quantum perturbation theory provides a faster way to derive the equations 
describing classical plasmas, once the ^ — )■ limit is taken (see [3] and references therein for 
a concise summary of the motivations leading to the quantum-mechanical approach). 

In recent years, special attention has been given to collective phenomena in quantum 
plasma systems [IHE], in which high densities of the components and small average inter- 
particle distances (of the order of the de Broglie wavelength) are assumed. In these works, 
a second quantization formalism is favored, leading to Wigner-Poisson and Wigner-Maxwell 
models. In references [H E], the authors compute the dielectric response function and elec- 
tronic conductivity of quantum and classical plasmas in the collisional regime using the 
Wigner-Boltzmann formalism [6], while in reference [9], the authors use the Schrodinger- 
Boltzmann formalism to write the dielectric function of a non- degenerate collisional plasma. 
In both formalisms, the collective response of the system is obtained through a perturbation 
on the chemical potential appearing in the usual Fermi-Dirac distribution for the equilib- 
rium electrons, and good agreement is shown between the classical results and those of the 
quantum models with /i — ;■ 0. 

Another field of application for the quantum approach is that of laser-plasma systems 
[TUHT^ . In particular, photon-plasma interactions have been discussed by many authors (see 
[3] and references therein). In our work, we compute the dielectric response function of a 
collisionless, one-component plasma system with neutralizing background, interacting with 
an external radiation field. To do so, we avoid second quantization arguments and use the 
Schrodinger description to account for the state of the plasma constituents. We extend the 
work developed in pTU], where there is no photon interaction contributing to the collective 
modes, in what they call a weak electron-radiation coupling; and that developed in [TT| WI\. 
where two limiting cases are discussed, i.e., the strong-field limit, in which only multiphoton 
processes are significant, and the weak- field regime where only single-photon processes are 
significant. 

For this study, the laser beam is treated as a classical plane electromagnetic wave in the 
dipole approximation. We consider the laser linearly polarized along the 2;-axes, with the 



electric field along the x-axes, taking into account a finite number of photons interacting 
with the electron plasma. This simple extension gives rise to a different dispersion relation 
for the electron waves. We see that an asymptotic value for the plasmon frequency exists as 
we increase the radiation wave number, and that this frequency never goes to zero, given a 
non-zero natural plasma frequency. 

This paper is organized as follows: in section II, we compute the state of the electrons 
through a unitary transformation [13] using the Schrodinger formalism; in section III, know- 
ing the fluctuations in the wave function due to the external potential, we calculate the 
dielectric response function; in section IV, we present the numerical scheme for obtaining 
the zeros of the dielectric function and discuss the results for the collective modes; in sec- 
tion V, we calculate the electric conductivity; in section VI, with an approximation on the 
plasmon frequency, we obtain an expression for the Landau damping term of the collective 
modes; we close with a summary of the main findings and difficulties. 

II. ELECTRON STATES 

For an electron under the presence of an electromagnetic wave, the Schrodinger equation 
takes the form 

i/*.(r.t) = m^. (1) 



with the Hamiltonian operator 



"=2;i^(P-'^W)^- (^) 



where p is the momentum of the electron, A(t) is the vector potential of the radiation 

A(t) = ( — sin(a;t) J x, (3) 

and uj is the frequency of the external radiation. To solve this equation, given a time- 
dependent potential, we use a unitary transformation [131 [H] of the form 



^(r,t) = [/<l>(r,t), (4) 

$(r,t) being the solution of the Schrodinger equation for a free electron and U the unitary 
operator given by 

f^ = exp ( ^a(t) -rWxpi ^I3{t) ■ p j exp ( ^r]{t) j , (5) 



where the functions (y.{t),j3{t) produce, respectively, translations in momentum and space, 
and rfit) is a phase factor. Substituting this into equation ([I| we obtain 

dihirA) i f da d/3 d?7\ ,, , ,,(9$(r,t) _, 



dt h\dt dt dt y ' ' dt 

Multiplying this equation by U\ we get 

*°=^ + i("W-eAW)V£-H<)-eAW) 

da d/3 da ^ d/3 d?7 ,„, 

where i^o is the Hamiltonian for the free electron. We proceed to solve the equation for a, /3 
and rj (linear terms in r and p and independent terms are set equal to zero) and find 

ipir, t) = exp ( TF{t) j exp (i7oA;a- (1 - cos(a;t))) exp(ik ■ r) exp ( -T^^t j , (8) 



for the wave function of an electron in an electromagnetic field given by (tsl), where we have 



simplified using F{t) = -2'jiujt + 71 sin(2u;t), 70 = -f^a , 7i = ^!^ and ek is the energy of 



the free electron with wave number k. 

We can see that (tsl) forms an orthonormal set. For a plasma, the presence of the external 
field generates local fluctuations in potential, and we may use (tsl) as basis to expand the 
wave function of these electrons in a local potential 

^,(r,t) = 5^ak(t)V'(r,t). (9) 

k 

Assuming a weak local potential of the form 

ip{v,t) = / dq / dr2exp(iq ■ r) exp(— ir2t)(y9(q, fi) + c.c. (10) 

we determine the coefficients a^it) using usual perturbation theory [T5] 

ak+q(^) = e ■ exp(-i7og:c) ^ i™Jm(ga;7o)v5(q, ^) 

exp (^ (^k+q — £k — ^^ ~ mhu) t) 



£k+q — ^k — ^^ — mhu! — iC 

where we use the identity 



(C^o+), (11) 



exp(iacosa;) = >^i™Jm(a) exp(ima;), (12) 



with Jm(«) being the Bessel function of order m. 
Finally, we can write the wave function as 

^fc(r, t) = es-^(*) exp (i7o/c^ (1 - cos(wt))) 

g-i^ki+ik-r J 1 + e ^ exp (-i7ogx cos(a;t)) 



Equation (13) is the wave function for the electrons of the plasma under the incidence of 



the radiation given by (pi). 



III. DIELECTRIC FUNCTION 



Knowing the states of the electrons via (13), we can obtain the fluctuation of the charge 
density 

p,(r,t) = -e|vl/,(r,t)|^-pl°)(r,t), (14) 

Pl. (r,t) being the charge density in the absence of a weak local potential, z.e, the charge 
density given by the ipi^^t) distribution. Neglecting terms of higher orders in ip, we have 

Pfc(r,t) = -e^ y^ e^'^'Xq, n) exp(-i7og^ cos(cut)) 

y'i°'Jm(gx7o) <^ TT^ T 

■^ I ^k+a — e-k — Ml — mnuj — iC 

m VIM 

gi{Q+mw)i ^ 

£^k+q - ^k + ^^ + m/icj + iC J 

Assuming a Maxwellian distribution /k for the electrons ( for a discussion on the choice 
of this distribution see ref [IS]), we have the total fluctuation as 

p(r,t) = ^/kPfc(r,t). (16) 



Using, once more, relation (12), we obtain 



p(r, t) = -e^ Yl e'''""'''V(q, ^) E i"~"' Jm(g.7o) Jm'(?x7o)e-'("^-'°')"*n(q, Q + mco), (17) 

q,f2 m,in' 



where 11 (q, Q,) = J2\i i — -iP-Mi-ic ^^ ^^^ electronic polarizability and m corresponds to the 



number of photons involved in the process [T7]. We take the real part of (17) to calculate 
the fluctuation. This fluctuation induces a potential in the medium given by the Poisson 
equation 

VVind(r,t) = -47rp(r,t). (18) 



Using (17) and the Fourier transform of (18), we obtain 



V5ind(q, fi) = ^^V5(q, fi) J^ J^(ga;7o)n(q, Vt + mw), (19) 

m 

which is the induced part of the full local potential 

^(q, n) = ^ext(q, ^) + ^ind(q, ^) = ^^^^^^, (20) 

where e(q, ^2) and (^ext(q, ^) are the dielectric function of the plasma and the external 
potential, respectively. The roots of the real part of this dielectric function give us the 
frequency of the longitudinal waves (collective oscillations, z.e, plasmons) in the plasma. 
We, therefore, separate e(q, Vt) into a real (cr) and an imaginary (ej) part. We proceed to 
the classical limit by letting /i — )■ and, after some algebraic work, we are left with 

e„(q,n) = l~.i^MM.-<.)l, • (l + ^ + ^^) e(-*). (21) 

III ^ / 

where A = mw + fi, ujp is the plasma natural frequency, e^ = 2'-fiu is the energy of the 
electromagnetic radiation and the average (f7(q, v)) is taken with respect to the function fk- 



IV. COLLECTIVE MODES 



The roots of (21) give us the frequencies of collective oscillations in the plasma as a 
function of their wave number (f2(q)). Here we encounter a numerical difficulty, as the 
sumation over m does not allow us to get an analytical result for this frequencies. However, 
if we fix a value for q, we obtain an expression depending only on Q, which is easily solvable 
using a bisection method in FORTRAN language [18]. The sum over m is truncated for 
a value of m beyond which the terms become neglegible. In this way, we take values of q 
from zero to about 20000m~^, and for each q we get a frequency value. So, we can plot the 
dispersion relation {Q versus q) of this system for various values of the radiation frequency. 



In the same manner, we can fix values of E and get tlie dependece of Q with E. These plots 
are shown in figures ([I| and ([2]). We see that, for any given value of to (the frequency of 
the radiation), the asymptotic value for the plasmon frequency is the same, the difference is 
that for larger values of u the plasmon frequency decays slower. From this we see that, for 
a radiation with large enough energy, the plasma remains unperturbed. 
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FIG. 1. (left) The natural plasma frequency for this case is 6.10^^s^^ and E = WV/m. Notice that, 
as we increase the radiation frequency, the dispersion decays slower, (right) There is a assymptotic 
value for the plasmon frequency as we increase both m an q, which seems to coincide with the 
natural frequency. The curve for m = reproduces exactly that of ref. [lOj . 



We can also vary the number of photons involved in the process and see how the curves 
react for a natural plasma frequency of 6.10^^s^^. This plots are shown in figures ([l]) and 

0- 

In figure M wee note a clear restriction to the range of frequencies allowed to the plasmons 

as m increases. So, as more photons participate in the process, more energetic the plasmons 

are (given the same value of q). For all plots, the temperature of the plasma is ksT = 

1,6.10-iV. 

The dispersion relation obtained for the free plasma subjected to a radiation field is 



eR(q,fi) = 



(22) 



1--JEJ 
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l + 5^,e 



\2 



(23) 
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FIG. 2. (left) The natural plasma frequency for this case is, also, G.lO^^s ^ and the value of q is 
generated randomly in the interval of figure (IT]). For the E dependence, we see an almost exponential- 



like decay, which agrees with the exponential term that appears in ( 21 ). (right) We see the tendency 



of the curve to become an exponential-like decay as we increase the number of photons in the 
process. 



This equation is quadratic in fi. We can easily check (23) to confirm that there are pairs 



of solutions of the form ±Q. Furthermore, this solutions must also be even functions of q. 

Using typical parameters of discharge plasma we plot in figure (tsl) the dispersion relation 
for a radiation field in the radiofrequency range. 

As expected, for lower plasma frequencies, the plasmon frequency assumes smaller values 
as well, retaining the assymptotic value at the plasma frequency. We must note that for 
large |g|, the code presents strong numerical fluctuations. Therefore the behavior of the 
curve in this region can not be attributed to any physical cause. 

If we assume smaller values for the plasma frequency, we can observe numerical instabil- 
ities for much smaller values of the plasmon wave number, as shown in the graphs of figure 

©• 

Taking only the first quadrant, we examine the dependence of the dispersion relation 
with the amplitude of the external field (figure ([s]) ). For large values of E, the curve gets 
rougher and the numerical instabilities happen at smaller values of q. Actually, it is not yet 
clear if this instabilities are due to numerical fiuctuations or to a break down of our linear 
assumption. 
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FIG. 3. Plot obtained with E = lOV/m and co = 3 ■ lO^s^^. Dispersion relation for various values 
of the plasma frequency. The modes excited in plasmas with lower natural frequencies have lower 
energies. 

To check the validity of our code, we can easily plot the value of eR{q,Q) to see if it is 
actually close to zero. We observe in figure (|6| that our code shows consistent results. For 
larger values of |g| and \E\, the results get further away from the actual roots of e/j(g, Q). 

We can, also, estimate the imaginary part of the dieletric function in the classical limit. 

3/2 ,,,2 



ei(q,fi) 



TT / rUe 

2 \k^ 



^e-^/'^-^^J^,(g7o)A.e- 



me(>^in/q) 



q 



(24) 



This expression will be usefuU in calculating the dynamical conductivity, in the spirit of 
in the next section. 



V. CONDUCTIVITY 



The conductivity of the system can be obtained from the expression 

Airi 
e{q,n) = l + —a{q,n) 

If we write the conductivity as a = ur + icri, we arrive at the relations 

aR(g,fi) = —ei{q,n) 



a,iq,n) = ^[l-en{q,n)] 



(25) 



(26) 
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FIG. 4. (left) Dispersion relation for ujp = l.ells ^ . (right) Dispersion relation for cjp = l.elOs ^ . 
(bottom) Dispersion relation for Wp = l.el9s^^. 
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FIG. 5. Dispersion relation for various external field amplitudes. 

The real part of the conductivity presents a strong q dependent peak, shown in figure [7] 
For a fixed ^2, we can see, in figure (IS]) an oscillatory behavior in the imaginary part of 
the conductivity for small values of E. A strong, q dependent, peak is followed by a damped 
region. We see that a\ goes to zero in the region of collective behavior. (See figure (|2])). 



As expected from eq. (27), a small change [20] in Vt does not change significantly the 
value of a I as can be seen from figure (|9]). 
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FIG. 6. Real part of the dieletric function for the coUective modes: (left) Fixed E = lOV/m, 
(right) Fixed q = 210m~^. 
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FIG. 7. Real part of the conductivity for various values of q. 



VI. LANDAU TERM 



So far we have assume a real value for the plasmon frequency Q. Nevertheless, the 
dielectric function was defined in a such a way that it can assume complex values (e = eR+iei) 
for real q and Q. For complex Q, e(g, Q) is to be interpreted as the analytic continuation of 
e{q,Q) from the real Q axis. [3]. 

We write the solution of the equation e = in the form 



Q = r^R + 17. 
11 



(27) 
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FIG. 8. Imaginary part of the conductivity for various values of q. 



If we assume 




-1 ' 1 ' 1 ' 1 ' 1 ' r- 



1.0x10"n 





100 


200 


300 


400 500 






r 


\ 




fi = 2.e11 
q = 6000 


6.0x10"- 




1 


\ 






00- 


J 


f 






- 



100 200 300 400 500 



E(V/m) 



FIG. 9. 0"! for slightly diferent values of Q. 



I7I « |1]r|, |ei| << IcrI 



(28) 



we can expand e = and neglect terms of higher order in small terms to obtain 



eR(g, Or) + i7 



an 



+ iei(g,fiR) = 0. 



(29) 



r2=r2 



R 
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Thus, we get two equations for the real and imaginary parts 

e(g, fipt) =0 ^ Used in previous sections 

7--^ (30) 

\ an )nn 

Using our results, and neglecting higher orders in (1/Am), we obtain 

for the imaginary part of the plasmon frequency, responsable for the damping (or growth) 
of the collective modes. 



VII. CONCLUSION 

The number os photons involved in the interaction between the plasma and the electro- 
magnetic radiation is of extreme importance. We see that a single photon restricts the range 
of frequencies allowed to the plasmons (see figure ([I])), and, as we increase the number of 
photons, an assymptotic value appears for these frequencies. 

The plasmon frequency decays very rapidly with increasing field amplitude (E) in an 
exponential-like curve (see figure(p])). As we increase the number of photons in the process, 
we no longer see an oscillation in the curve of Q versus E as reported by [ID] . 

We successfully derived an equation for the conductivity of the system and a first ap- 
proximation for the imaginary part of the plasmon frequency. 

Although the codes present strong numerical fluctuations in regions of high valued electric 
field amplitude, the divergences of the dispersion relation in these regions can, also, be 
associated with the break down of the linear approximation (as well as the non-relativistic 
one). The effects of nonlinearty are under study and we are currently extending these results 
to compute magnetic field effects. 
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